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Order parameter fluctuations for the two dimensional Ising model in the region of the critical 
temperature are presented. A locus of temperatures T*(L) and of magnetic fields _B*(L) are identi- 
fied, for which the probability density function is similar to that for the 2D-XY model in the spin 
wave approximation. The characteristics of the fluctuations along these points are largely indepen- 
dent of universality class. We show that the largest range of fluctuations relative to the variance 
of the distribution occurs along these loci of points, rather than at the critical temperature itself 
and we discuss this observation in terms of intermittency. Our motivation is the identification of 
a generic form for fluctuations in correlated systems in accordance with recent experimental and 
numerical observations. We conclude that a universality class dependent form for the fluctuations 
is a particularity of critical phenomena related to the change in symmetry at a phase transition. 

PACS : 64.60.Cn, 05.40.-a, 05.50.-Fq 



I 

o 
o 



> 

i> 
in 
in 
o 

o 



X3 
O 

o 



X 



INTRODUCTION 

There is ubiquitous interest in systems with extended 
spatial and temporal correlations, from all areas of 
physics. Recently, observations of fluctuations in global, 
or spatially averaged measures of such correlations have 
provided a possible link between critical phenomena 
and non-equilibrium systems. That is, the probability 
density function (PDF) for order parameter fluctuations 
in the low temperature critical phase of the 2D-XY 
model, shown in Figures 1 and 2 below, is very similar 
to PDFs 0] for fluctuating spatially, or temporally 

for 
for nu- 
and of 




averaged quantities in turbulent flow 0, 
electro-convection in nematic liquid crystals_ 
merical models of dissipative systems |y, 0, 
"self-organized criticality" fluctuations of river 

heights [ij [l^ , as well as for other equilibrium systems 
close to criticality [H E [11 [H. The simplicity of 
the 2D-XY model allows a complete understanding of 
fluctuation phenomena in this case [lil IitI H^ . The 
contrary is true for non-equilibrium systems; the lack of 
microscopic theory makes the problems extremely com- 
plex. Phenomenological observations and analogy with 
better understood systems can therefore be extremely 
useful. 

However, for other critical systems, such generic be- 
haviour seems only to be observed under restricted 
conditions. For example, it has been shown that mag- 
netic fluctuations in the two-dimensional Ising model, 
at a temperature, r*(L) (where L is the system size 
in units of the lattice constant), below but near the 
critical temperature T^., are similar to those of the 
2D-XY model US 111 111 El- Given the important 
role played by universality classes in critical phenomena, 
this is rather surprising. In fact, it is well established 
that critical fluctuations, as measured at Tc [l9l l2oll2l| . 



depend, in general on the universality class of the model, 
as well as on the shage ^1] of the sample and on the 
boundary conditions (22l | . The apparent similarity of the 
form of the fluctuations, over and above the universality 
class of the models under consideration therefore seems 
rather puzzling (13j . 

Given the generality of the above observations in more 
complex systems, it is important to understand this 
point. In this paper we address it in detail for the 
two-dimensional Ising model and in doing so pose the 
following questions: Is the similarity in form quantita- 
tive, or only qualitative? Is this PDF really a measure 
of critical fluctuations in the problem? Finally can we, 
from this investigation shed any further light on the 
reason for the apparent " super- universality" observed in 
the wide range of experimental and numerical systems ? 
In answer to these questions, we show that the distri- 
bution functions for the Ising and XY models at r*(L) 
are similar, with the latter representing an excellent 
fit over almost any accessible window of measurements 
for experimental systems. However they are not the 
same functions. Numerical evidence suggests that the 
difference will remain in the thermodynamic limit. The 
origin of the difference is the structure of phase space 
associated with the Ising transition. For large amplitude 
fluctuations the system, localized in one half of phase 
space is able to surmount the barrier separating it from 
the other, symmetric half of phase space. If instead 
one approaches the critical point by applying a small 
magnetic field, symmetry remains broken and one finds 
a field B*{L) giving excellent quantitative agreement 
between the PDF for the Ising model and the 2D-XY 
model. 

We confirm that, despite the small value of the correla- 
tion length, ^ [ij, [l^ at T*(L), the system does show 
evidence of correlations on all scales up to a length of 
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the order of the system size. A consequence of this is 
the development of coherent structures, the clusters of 
spins, up to this macroscopic scale, that dominate the 
exponential tail of the distribution. 

From the above, we conclude that the dependence of the 
PDF on the universality class comes from the structure 
of phase space and the way in which symmetry is 
restored on passing through the phase transition. If long 
range correlations develop, without the insuing fluctua- 
tions exposing a structured phase space, such as occurs 
in the Ising transition, then the form of the fluctutations 
will be largely independent of universality class. This is 
the case in the 2D-XY model and wc propose that it 
is key to the approximate " super- universality" observed 
in a large array of correlated systems. 
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ORDER PARAMETER FLUCTUATIONS IN THE 
2D ISING MODEL 



The Hamiltonian for the Ising model is given by: 



J ^ ^ SiSj^ Si — il 



<hJ> 



with J > 0, the exchange constant. Wc study the model 
on a square lattice of size TV =L x L , with periodic bound- 
ary conditions and with lattice parameter a ~ I. The 
order parameter is defined as the modulus of the magne- 
tization: 
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Previous studies, by Bramwell et al. [lo| . by Zheng and 
Trimper and Zheng suggest that there is a tem- 
perature, T*(L) just below the critical temperature, for 
which the PDF is close to that of the 2D-XY model in 
the low temperature phase. It was further suggested that 
T*(L) scales towards the critical temperature as the ther- 
modynamic limit is taken and can thus be interpreted as 
a critical phenomenon 0|. The first step in our study 
was to confirm this result. Note that, within the spin 
wave approximation, the PDF for the 2D-XY model has 
the form shown in Figures ^ and |21 independently of 
temperature and hence of critical exponent jl6l |. The 
PDF was calculated using the Swendscn-Wang Monte 
Carlo algorithm [2^ for various sizes L , between 32 
and 512. Here we define T*(L) as the temperature for 
which the skewness of the distribution, 7(T*(L)) is equal 
to that for the 2D-XY model, 7xy, for data restricted 
to the window /i = (m — {m))/a G [—6; 3]. We are 
able to establish this criterion with a numerical preci- 
sion 7(r*(L)) = 7xY ± 0.01, giving an error for T*(L) 
of less than 2%. This criterion, though precise, is arbi- 
trary and we could choose others. However, as it is only 



FIG. 1: Order parameter PDF for T = 2.33J = T^L = 64) 
win, r*(L = 64) = 2.11J (x) and T = 1.54J (o) , with 
L — 64. The inset is a zoom on /i G [—5; 3]. 




FIG. 2: PDF for a subsystem L =32, in system of size Lo = 
128, for B = B*{1 = 32) = 0.0035J at T^{1 = 32). The plain 
curve is the PDF for the 2D-XY model. 



an approximate agreement between the two sets of data, 
making alternative definitions does not alter the results 
discussed below. Figure shows \og{<Tll{m)) against fi 
for various temperatures and for L =64. The solid line is 
the PDF for the 2D-XY model The best fit is for 

r*(L = 64) = 2.11J, while the critical temperature for 
the infinite system, as given by Onsager's exact solution 
is Tc = 2.25 J. One can see that r*(L) is significantly 
shifted below Tc- At higher temperatures, the probabil- 
ity for fluctuations below the mean is a concave function 
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of and there is no region approximating to an expo- 
nential tail. The PDF is cut off at a finite value of /i 
corresponding to the constraint m = 0. As one expects, 
this corresponds to a turning point in probability and re- 
flects the access of the finite size system to the complete 
and symmetric phase space. 

A distribution with an exponential or quasi-exponential 
tail cannot correctly describe this minimum and there- 
fore for symmetry reasons the distribution for the 2D- 
XY model cannot exactly describe the data for the Ising 
model in zero field, for any temperature below T'c(L) |24| . 
However, at T*(L) the fit is good for fluctuations out to 
—5(7 below the mean magnetization, corresponding to a 
probability density of 10~^ |1!tI |. From an experimental 
point of view, reliable data for fi < —6 would be excep- 
tional |3i la Id I and in this sense the agreement between 
the Ising and the XY model data is very satisfactory (see 
inset Figure 1). Here we have good statistics for fluctu- 
ations out as far as —8a from the mean, from which an 
upturn away from the exponential tail is evident. This is 
a consequence of the extremum in PDF at m = 0. The 
effect is independent of system size for the values stud- 
ied and within the numerical error obtained, although 
we cannot exclude the possibility of corrections to scal- 
ing that disappear slowly on the scale of the sizes stud- 
ied. For temperatures below r*(L) the turning point at 
771 = is moved outside the accessible window of mea- 
surement, but the PDF is not sufficiently skewed to give 
a good fit to the XY model data. For temperatures be- 
tween T*(L) and T = 1.54,7 the the statistics are poor 
for large fluctuations, which is consistent with having just 
a few rare events taking the magnetization out towards 
the constraint m ~ 0. At lower temperature still, the dis- 
tribution crosses over to a Gaussian, as expected for an 
uncorrclatcd system. The relevance of the turning point 
in the PDF in the critical region depends on the univer- 
sality class and it is one of the ways in which universality 
class dependent critical fluctuations appear. 



BOUNDARY CONDITIONS AND MAGNETIC 
FIELD 

It is clear that the quality of the comparison with the 
2D-XY model would be improved if the turning point in 
the probability density at m = was either displaced 
or removed. Two obvious ways in which one might do 
this are firstly in changing the boundary condition and 
secondly in adding a magnetic field. Changing from pe- 
riodic to fixed, or window boundaries one can expect to 
observe small changes in the form of the universal scaling 
function [lil 12^ 1 . One might think that these boundaries 
would make a finite size system more rigid with respect 
to a global spin flip, thus reducing the probability of a 
microstate with m = and improving the quality of the 
fit. We have studied the distribution for a window of size 



L embedded in a larger system of size Lq . The fits to the 
PDF for the 2D-XY model are qualitatively better, but 
the same upturn away from the exponential tail is ob- 
served for fluctuations of more than about 5a from the 
mean. The situation is not quantitatively changed com- 
pared with periodic boundaries. The data are not show 
here. 

A real quantitative improvement is found however, for 
the study of fluctuations in a magnetic field. Approach- 
ing the critical point by fixing T = Tc(L) and applying 
a small field B, a field B*{L) is defined for which the 
PDF of magnetic fluctuations gives the best fit to the 
2D-XY data. The data for B*{L) = 0.0035 J for win- 
dow boundaries, with L = 64 and Lg = 128, are shown 
in Figure |2| To the eye the quality of the fit is excel- 
lent. This is confirmed quantitatively by measuring the 
skewness, 7, and kurtosis k, of the distribution. We find 
7 = 0.890±0.01, K = 4.495±0.01, which are in excellent 
agreement with those of the 2D-XY model (7xy = 0.890 
and KxY = 4.489). The application of a magnetic field 
breaks the symmetry and removes the minimum from the 
PDF for the order parameter. This can be seen in Figure 
where we show the evolution of the PDF, for a win- 
dow of size L = 32 in a larger system of size Lq = 128 for 
field varying between zero and B*{L ^ 32). For zero field 
the minimum in probability at to = is clearly visible 
and the fluctuations below the mean are cut off by the 
constraint to ^ 0. As the field increases the mean mag- 
netization increases and the variance, a reduces, pushing 
the constraint to ^ out to larger negative values of 
fi. At the same time the end point of the distribution 
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FIG. 3: PDF for a subsystem L =32, in system of size Lo = 128, 
for various values of the magnetic field B at Tc(\- = 128). The 
plain curve is the PDF for the 2D-XY model. 
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is no longer a minimum, that is, the PDF terminates 
with a finite slope. Within the window of observation, 
one clearly sees an exponential tail developing which ap- 
proaches asymptotically that of the 2D-XY model. For 
larger fields the asymmetry reduces and the PDF crosses 
over to a Gaussian. The curve for the 2D-XY model 
therefore seems to define an envelope giving the maxi- 
mum possible asymmetry as the field varies in the region 
of the critical point. Behaviour for periodic boundary 
conditions is similar. One again sees the development of 
an exponential tail with the same slope as the 2D-XY 
model, but for the best fit the cut off corresponding to 
TO = remains within the window — 8 < /i < 3. The data 
are not shown here. 



CORRELATION LENGTH AT r*(L) 

In this section we concentrate on the critical properties 
of T*(L). Figure^ shows how T*(L) scales with system 
size. As noted in references [Ulia, r*(L) scales with 
L as : Tc - T*(L) ot L"^/", with i/ ~ 1, the expected 
value conforming to the scaling hypothesis for the 2d 
Ising model. These finite size scaling results imply that 



The numerical data is fitted well by the expression 
GL(r,T)-i-e-^, 

with rj = 0.24±0.01, in good agreement with the theoret- 
ical value rj = J. Values of ^(L) and ^(L)/L arc shown in 
tabled The main criticism of Zheng and Trimpcr |l3l| . 
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FIG. 5: Correlation functions for L = 128 and various tem- 
peratures T = 2.30 J(+), 2.27J(x), 2.24J(*), 2.17J(n). The 
plain curves are the best fits. 




FIG. 4: Evolution of T*(L) with L . The solid line is the best 
linear fit. 

the magnetic correlation length ^(L) is diverging with 
the system size at T'*(L), such that C('-)/L = constant. 
This is confirmed by the observed universality of 11 (/x) 
along this locus of points The correlation length 

can be computated from the spin-spin connected corre- 
lation function: 

GL(|r,-r,|) = {S{r,).S{r,))-{S{r,)).{S{r,)). 

Different curves, obtained for L =128 and various tem- 
peratures are plotted on Figure [S] 



that ^(L) is small, is immediately apparent, however the 
ratio ^(L)/L ~ 0.03 is indeed constant, to a good approx- 
imation, confirming that the correlation length is diverg- 
ing with system size, as conjectured. In this respect we 
are indeed dealing with a critical phenomena despite the 
small values measured for the system sizes considered. 
This "small, yet diverging" quantity plays an important 
role in the approximate universality observed between 
disparate systems. 

CLUSTER SIZES AND DISTRIBUTION 

Definitions of clusters in the Swcndscn-Wang algo- 
rithm is a good way to study structures in the 2d Ising 
model A cluster is defined as a connected graph of 
spins in the sense of reference j2fij . with magnetization 
opposed to the spatially average value over all the lat- 
tice. In concrete terms, a cluster is a white object in the 
snapshots shown in Figures El [7| and |H1 Cluster sizes are 
calculated for each generated spin configuration and av- 
eraged over many realizations. 



L 


32 


36 


40 


44 


52 


56 


64 


128 


m 

C(L)/L 


0.83 
0.026 


1.0 
0.028 


1.1 
0.027 


1.25 
0.028 


1.45 
0.028 


1.42 
0.025 


1.7 

0.026 


2.9 
0.026 



TABLE 1: i and ^/L for various L . 



Given the small value for the correlation length at T*(L) 
one might expect the range of cluster sizes to be ex- 
tremely limited. One of the surprises of this study is 
that this is not the case. One can get a feeling for this by 
first studying snap shots. In Figures IHl [3 and |H1 wc show 
three configurations, the first has magnetization m close 
to the mean value. One can clearly see a range of cluster 
sizes up to a characteristic size that is small compared 
with L. The second snap shot shows a configuration with 
magnetization four standard deviations below the mean, 
m = (m) — 4cr. A much bigger cluster is present. The 
large fiuctuation is due to the presence of this large coher- 
ent structure. This is very different from what one would 
expect for fluctuations of an uncorrelatcd system: in this 
case, a large deviation would correspond to a configura- 
tion with many, small and uncorrelated clusters appear- 
ing spontaneously. This scenario is extremely unlikely, 
which is why fiuctuations away from the critical point 
are Gaussian. Here it is clear that, while the correla- 
tion length fixes the size of typical clusters, much bigger 
clusters are not excluded. They are rare events, but not 
so rare as to be experimentally irrelevant. This can be 
compared with the configuration taken at T'c(L), shown in 
Figure IHl Here, coherent structures spanning the entire 
system are not rare. In this situation the PDF depends 
strongly on the universality class and generic behaviour 
is not expected. Moving along the locus of points r*(L), 
the size of typical clusters scales with L through the scal- 
ing of ^(L) and we expect the size of rare clusters to scale 
in the same way. In this sense the only length scale in the 
problem along T*(L) is L. Apart from this, the system 
is scale free, despite the small values of ^(L) extracted 
numerically. 

From the mapping of the Ising model at criticality onto 
a percolation problem l27l | one expects the clusters 
to have a power law distribution of sizes. Results for the 
distribution of clusters sizes P{s), obtained for L = 128 
for various temperatures below Tc{\-) are presented in 
Figure]^ As the temperature is below Tc{\-) the system 
possesses a spanning cluster whose spin direction defines 
the global magnetization direction. Note that, while for 
standard Metropolis Monte Carlo the magnetization di- 
rection of the spanning cluster changes extremely rarely, 
for Swendson-Wang Monte Carlo symmetry is not broken 
and the direction oscillates between "spin up" and "spin 
down" . The spanning cluster is included in the statis- 
tics shown in Figure El As expected, as the temperature 
approaches T'c(L) the probability of finding large clusters 
increases and the distribution approaches a power law. 
However, for temperatures below r*(L) there is a sepa- 
ration of scales, with a gap in probability between the 
largest non-spanning cluster and the spanning cluster. 
The gap closes at about T*(L) and above this temper- 
ature the statistics of the spanning and the secondary 
clusters are mixed 0| . When this is the case the largest 
clusters of up spins and down spins will be of the same 




FIG. 6: Typical configuration at r*(L). 





FIG. 8: Typical configuration at Tc. 
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FIG. 9: Clusters size distribution P(s) for L =128 and T = 
2.29J(+), r*(L = 128) = 2.25J(x), T = 2.22J(*) and T = 

2.13J(n) {S = iVclustor/iV). 




FIG. 10: Clusters size distribution at r*(L) for L =64(A), 
128(o), 256(+) and 512(x). The solid lines are the best power 
laws obtained. 



size and there will be a non-negligible probability of hav- 
ing zero magnetization. The mixing of statistics of the 
spanning and secondary clusters is therefore perfectly 
consistent with the observation that the turning point 
of the PDF. at m = 0; moves into the window of numer- 
ical or experimental observations for T above r*(L) and 
that the tail of the distribution n(/i) is no longer well 
approximated by an exponential. 

At T*(L) the cluster distribution is well represented by 
a power law out to cluster size s — ^^'j^^"*" < 



m- This 



is just about the size of the large cluster in the second 
snap shot of Figure [7| Above this size, corrections to 
scaling are manifest as one might expect, but it is worth 
noting that the limit of power law behaviour corresponds 
to clusters of sizes well in excess of the correlation length 
extracted from the autocorrelation function. As shown 
in Figure [TUl the range of application of a power law dis- 
tribution increases with system size. For L — 512 it ex- 
tends over almost six orders of magnitude of probability 
and three of cluster size. The fitted exponent decreases 
with L, as shown in table UTI For L = 512 the best fit 
is for P(s) ~ s-^C-) with r(L) = 2.2 ± 0.1. This seems 
to correspond to the estimated exponent for percolating 
clusters U^E^,t = 2 + ^c^ 2.1, where D/ = 187/96 
is the fractal dimension for the 2D Ising model [s^l , 



L 


16 


32 


64 


128 


256 


512 


r(L) 


3.6 


3.6 


2.5 


2.4 


2.3 


2.2 



TABLE II: Exponant r(L) of the power laws obtained for the 
cluster size distributions. 



confirmation of the critical nature of our observations. 
In this section we have shown that clusters at T*(L) 
show, rather remarkably, both scale free behaviour, and 
a separation of scales: Typical cluster size, maximum 
secondary cluster size and spanning cluster size are all 
fixed by L, however their amplitudes are sufficiently dif- 
ferent to ensure a minimum of interference between the 
three scales. In the low temperature phase of the 2D-XY 
model, the ratio of variance to mean magnetization is, 
cr/ (to) = AT /J [13, where A « 0.04. This is a crit- 
ical phenomenon, in that the ratio is system size inde- 
pendent. However, it goes to zero at zero temperature, 
meaning that the critical fluctuations have zero ampli- 
tude and despite their singular nature, the system visits 
an infinitesimal part of phase space near the mean value 
of the magnetization. The separation of scales for cluster 
formation in the Ising case is analogous to this. It corre- 
sponds to a situation where the fluctuations are critical, 
but where the limits of phase space (to = 1, m = 0, 
or to = —1) are not approached. We propose here that 
observed generic behaviour and universality classes inde- 
pendance must have its origin in this point. 



EXTREME CLUSTERS DISTRIBUTIONS 

The form of the distribution in Figure 2 bares a strong 
resemblance with Gumbel's first asymptotic solution for 
extreme values 0, 0| . Indeed there have been a series 
of recent papers searching for a connection between this 
form for global fluctuation and extreme statistics [i^ls^f . 
Such a connection has failed to emerge in Gaussian in- 
terface models, related to the 2D-XY model 0|, but 
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as there is here direct access to obvious real space ob- 
jects, the clusters, it seems natural to investigate extreme 
statistics for the clusters. One relevant question is: does 
the largest cluster dominate Il{fi) rendering the problem 
of non-Gaussian fluctuations an extreme value problem 
for the clusters ? The answer is no! While the distribu- 
tion for the largest cluster is skewed in the right direc- 
tion and looks qualitatively quite similar to Ii{fJ-), it is 
not sufficient to reproduce the global fluctuations quan- 
titatively. This conclusion has been tested in detail: We 
express the magnetization as 

oo 

N ' 

with Uj = 0,Vj > jniax- Approximate order parameters 
TOfc arc constructed: 

fc 

En.j 

If extreme values statistics are relevant for the complete 
order parameter PDF, then starting from k = I H^fij;), 
with /ifc — {mi; — (mi;)) /a should converge to n(/i) for 
just a few values of fc. Results arc shown in Figure [TTI for 
fc = 1,5, 10 and compared with the complete order pa- 
rameter PDF. The convergence is slow and fc = 10 is not 
sufficient to reproduce well the global PDF. We conclude 
that all the clusters arc required to reproduce the global 
fluctuations. This is not then an extreme value problem, 
at least in terms of clusters. 

This result is rather similar to that obtained for the 2D- 
XY model 0| and for related Gaussian interface mod- 
els although in the latter more detailed information 
on the microscopic distributions is available. The non- 
Gaussian fluctuations occur because the largest cluster 
makes a macroscopic contribution to the many body sum. 
However, this contribution docs not dominate, it is the 
same order as the sum over all the other clusters. In this 
sense critical fluctuations are a marginal case between 
statistics dominated by the majority, leading to the cen- 
tral limit, and statistics dominated by a single event, for 
example Levy statistics. 

INTERMITTENCY AND CRITICAL 
FLUCTUATIONS 

Figure^lshows a sequence of magnetization values for 
configurations generated by a Metropolis Monte Carlo al- 
gorithm at temperature r*(L). The asymmetric nature 
of the distribution is evident from the stochastic time se- 
ries. We remark that the data bare a striking resemblance 
to the time series of injected power into a closed turbu- 
lent flow at fixed Reynolds number Q . The two systems 
share the characteristic of making large deviations from 




FIG. 11; PDF of rrik, for k = l(-h), 5(o), 10(x), and non 
approximate order parameter (solid line), for L = 64 at 
r*(L = 64). 



the mean value, on a scale set by the variance of the dis- 
tribution a. This is the so called " intermittency" which 
is so important in studies of turbulence. It is not the 
purpose of this paper to discuss the detailed properties 
of intermittency, we simply make a comment concerning 
the scale of intermittency in the 2D Ising model in the 
region of the critical point. In Figure El ^ similar time 
series is shown for a simulation at Tc(L). One can clearly 
see that the range of fluctuations, on the scale of a is 
smaller at Tc(L) than at T*(L). This is an unconven- 
tional way of characterizing fluctuations in the context 
of critical phenomena. Usually one discusses the suscep- 
tibility, X ~ (l/^)'''^! which is a diverging quantity at 
a critical point. For a flnite size system, rc(L) is gen- 
erally defined at the temperature where x and hence a 
is a maximum [s^]. What we observe here is that, al- 
though the fluctuations on an absolute scale are maxi- 
mized at T'c(L), on a relative scale, as fixed by a they 
are maximized at T*(L). Between r*(L) and T'c(L), large 
deviations from the mean value are cut off because of 
the constraints on the phase space in which the fiuctua- 
tions can occur. Specifically, the free energy barrier be- 
tween "spin up" configurations and "spin down" config- 
urations is surmountable with the consequence that the 
limits < (m) < 1 become the relevant limits for fiuc- 
tuations and the variance and the mean order parameter 
values become of the same size. In this way the mag- 
netic symmetry, which is broken in the ordered phase, is 
re-established on passing through the phase transition. 
For temperatures up to T'*(L), despite the growing fiuc- 
tuations the system behaves almost as if only a single 
phase existed, for which the order parameter could vary 
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FIG. 12: Metropolis evolution of m at r*(L). 
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FIG. 13: Metropolis evolution of m at Tc{l). 

within the limits — oo < m < oo. The distribution for 
the 2D-XY model corresponds to the case where this is a 
hmitingly good approximation |23l | . If, therefore we make 
a rather loose definition of intermittency, as the tendency 
to make fluctuations away from the mean, on the scale of 
a with a probabihty that differs from that for a Gaussian 
fimction, then maximum intermittency occurs at r*(L) 
and not at Tc(L). 

CONCLUSION 

In this article we have compared the order parameter 
fluctuations of a finite size 2D Ising model, in the region 
of the critical temperature, with those of the 2D-XY in 
the low temperature phase. Our motivation for doing 
this is the observation that fluctuations of global quanti- 



ties in a wide range of different correlated systems are of 
similar form to the 2D-XY model, or equivalent Gaussian 
interface models jl^]. Approaching the critical temper- 
ature Tc, either along the zero field axis, or applying a 
small field, B, at the critical temperature one can iden- 
tify a temperature T*{L){B = 0) and a field B*{L), close 
to the critical point for which the order parameter fluctu- 
ations are similar to those of the 2D-XY model. We have 
established the critical scaling behaviour of the locus of 
temperatures T*(L). The correlation length ^ diverges 
with system size along the locus of temperatures r*(L) , 
as one might expect for critical scaling. However, the 
ratio ^ 0.03 is a small constant. This is a key point: 
Critical fluctuations are usually associated with a phase 
transition and a change in phase space symmetry. In the 
2D Ising model the symmetry between the two compet- 
ing phases is reflected in the form of the PDF, imposing 
that m = corresponds to a turning point in the prob- 
ability density. However the 2D-XY model is an excep- 
tion, in that there is a continuous line of critical points 
in the low temperature phase, but no phase transition or 
associated change in symmetry |35j . To an excellent ap- 
proximation the critical fluctuations occur in an uncon- 
strained phase space jl^. Hence, criticality in the Ising 
model can only resemble that in the 2D-XY model if the 
change in symmetry is not apparent. This corresponds 
to the condition that the ratio ^/L is small. In the case 
of B*{L), the agreement between the PDFs for order pa- 
rameter fluctuations in the two systems is exceptionally 
good: High quality numerical data for the Ising model 
arc indistinguishable from the analytical results for the 
2D-XY model. However, differences arc observable for 
data at T*(L). The agreement is better in the former 
case, as the field breaks the symmetry between the two 
competing ordered phases, thus eliminating the turning 
point, or the minimum value from the PDF. In this case, 
the phase space available for fluctuations strongly resem- 
bles that of the 2D-XY model. 

Driven, non-equilibrium systems, showing strong cor- 
relations, such as turbulent flow 0, 0|, resistance net- 
works 0, self-organized critical systems 0|, or grow- 
ing interfaces, resemble the 2D-XY model, in that there 
is no associated phase transition and no sudden change 
of limits that constrain the fluctuations as the correla- 
tions build up. We propose that observation of a non- 
Gaussian PDF with finite skewncss and an exponential 
tail, for fluctuations in a global quantity, is a characteris- 
tic of correlated fluctuations in an effectively unbounded 
phase space. The example of the 2D Ising model serves to 
show that this is not the case for critical phenomena asso- 
ciated with a 2"*^ order phase transition. However, this is 
a detail specific to critical phenomena, related to the fact 
that the fluctuations become so large that they allow the 
system to explore the whole allowed phase space of or- 
der parameter values. From this analysis one could argue 
that fluctuations at T*(L) are not strong. This is true on 
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an absolute scale: For example the susceptibility, which 
is a measure of the variance of the distribution, is small 
at r*(L) compared with the maximum value from which 
one defines Tc(L). However at this temperature extreme 
fluctuations compared with the variance are capped by 
the constraint, < to < 1. Parameterizing in terms of 
the reduced variable fi = (m — {m))/a changes this con- 
clusion: Fluctuations in /x are essentially unbounded at 
T*(L) while they are constrained at T'c(L). Hence the 
largest fluctuations in ^ occur at r*(L) and not at T'c(L). 
ft is the variance of the distribution which defines the 
scale of the fluctuations that one might observe exper- 
imentally or numerically, and so in this sense r*(L) 
is highly relevant. Physically, this means that at T*(L) 
large fluctuations arc rare enough not to modify a, but 
not too rare to be observed. As long as the constraints 
relevant to the phase transition remain unimportant, in- 
creasing the level of fluctuations, that is increasing the 
ratio (to) /a there will be little, if any evolution of the 
PDF. This is exactly what is observed for the 2D-XY 
model: The PDF has the universal form discussed above, 
independently of the critical exponent along the critical 
line 0, [i^. This seems consistent with the observed 
generic behaviour in disparate systems and also appears 
to be compatible with recent renormalization group cal- 
culations on Gaussian interface models with quenched 
disorder [s^ . In the latter a study is made of the PDF 
for Gaussian interfaces (of which the low temperature 
phase of the 2D-XY model is an example) in the pres- 
ence of disorder. This is shown to be highly irrelevant, 
with the PDF being unchanged from that of the under- 
lying Gaussian model. 

It remains to quantify what we mean by "little evolu- 
tion of the PDF" . In this paper we have shown that 
at r*(L), or B*{L), dependence on universality class is 
largely absent. It has been further shown that, while the 
boundary conditions are important for quantitative com- 
parison, there effects are not very significant. However, 
we have remained firmly in two-dimensions. Moving to 
three dimensions will no doubt lead to variations and this 
will prove an interesting test for our explanation of the 
observed approximate universality for global fluctuations 
in correlated systems. 
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